Generalized Kubo formula for spin transport: A theory of linear response to 

non-Abelian fields 



o 
o 

(N 

wo: 
< 

in 



CZ2 



13 

o 
o 



> 

\o 
o 
in 
o 

O 

-i— > 
c3 



i 

C 

O 

o 



X 



Pei-Qing Jin and You-Quan Li 

Zhejiang Institute of Modern Physics and Department of Physics, 
Zhejiang University, Hangzhou 310027, People's Republic of China 
(Received February 6, 2008) 

The traditional Kubo formula is generalized to describe the linear response with respect to non- 
Abelian fields. To fulfill the demand for studying spin transport, the SU(2) Kubo formula is derived 
by two conventional approaches with different gauge fixings. Those two approaches are shown to be 
equivalent where the nonconservation of the SU(2) current plays an essential role in guaranteeing 
the consistency. Some concrete examples relating spin Hall effect are considered. The dc spin con- 
ductivity in response to an SU(2) electric field vanishes in the system with parabolic unperturbed 
dispersion relation. By applying a time-dependent Rashba field, the spin conductivity can be mea- 
sured directly. Our formula is also applied to the high-dimensional representation for the interests 
of some important models, such as Luttinger model and bilayer spin Hall system. 

PACS numbers: 72.25.-b, 72.10.-d, 03.65.-w 



I. INTRODUCTION 

Kubo formula, one of the most important formulas in 
the linear response theory, has been widely used in con- 
densed matter physics since it was derived by Kubo [l| 
for the electrical conductivity in solids. There are several 
kinds of Kubo formulas for the external fields to which 
the system responses are different. However, these for- 
mulas, such as those for the electrical conductivity and 
the susceptibility, all describe linear responses to the U(l) 
external fields. 

Recently, a newly emerging field, spintronics has 
absorbed much attention for its promising applications in 
quantum information storage and proces sing . Sp in Hall 

effect 0, IE IE S H, 13, E3, Ell E3, 13 El, US EE (Hi , as 

a candidate method to injecting spin current into semi- 
conductors, is also discussed intensively. In this effect, 
the spin-orbit coupling is necessary no matter intrinsi- 
cally or extrinsically. As this coupling can be re gard ed 
as a contribution of the SU(2) gauge potential [l8|, a 
new version of linear response theory in SU(2) formu- 
lation is necessary. Nevertheless, most of the previous 
works have mainly focused on the linear response of such 
system to an external electric field, and hence the tradi- 
tional Kubo formula was adopted directly except Ref E3 
which dealt with non-Abelian response and considered 
spin Hall effect in the presence of an SU(2) gauge field. 
There were some papers [2(| discussing the responses 
to a spin-orbit coupling with spatially varying strength, 
but the authors employed other approaches rather than 
Kubo formula as the SU(2) Kubo formula has not been 
established. It thus becomes inevitable to develop a gen- 
eralized Kubo formula so that the linear response to the 
external SU(2) gauge fields can be evaluated. 

In present paper, we derive a formula which describes 
the linear response to an SU(2) external field using the 
strategy ever employed by Kubo for the U(l) case. It 
is not a straightforward derivation since the algebra is 
totally different. Especially, the expression of the SU(2) 



"electric field " evolves one more term of gauge poten- 
tials [12 than the U(l) case due to its non-Abelian fea- 
ture. It seems obscure to directly find the equivalence 
between the Kubo formulas derived with different gauge 
fixings. We will show that the extra term in the SU(2) 
"electric field " precisely corresponds to the nonvanishing 
term in the "continuity- like" equation E3 which includes 
the spin procession [2, HI- Its origin stems from the 



definition of the conserved current |18| in the presence 
of the SU(2) field. Since one of recent research interests 
focuses on the spintronics, some explicit examples in spin 
Hall systems are discussed in terms of our SU(2) Kubo 
formula, such as the spin susceptibility and current in 
response to the effective spin-orbit coupling [l3l Hij . In 
spin Hall effect, the spin conductivity is believed to be 
canceled by the effect of disorder in two-dimensional elec- 
tron gas |S|. It is due to the parabolic unperturbed dis- 
persion relation [2o|. In such a system, the spin current 
in response to the external spin-orbit coupling also van- 
ishes. Then the systems with nonparabolic unperturbed 
dispersion relation become significant. In such systems, 
the spin conductivity in response to either U(l) or SU(2) 
external fields does not vanish. An experimentally acces- 
sible case is also given in which the spin conductivity is 
related to the dielectric function. We also extend the ap- 
plication of our formula to a high-dimensional represen- 
tation, saying spin-3/2 representation, which is related to 
some important systems, such as the Luttinger model p3 | 
and the bilayer systems |2J| . The spin conductivity in the 
Luttinger model vanishes, which describes the response 
to the effective field of structural inversion asymmetry. 

The paper is organized as follows. In Sec.[Hl we derive 
a general Kubo formula with respect to a single- frequency 
SU(2) external field at zero temperature. Then we show 
in Sec. II I II that this formulation is consistent with the 
one by choosing a zero-frequency external field at the 
very beginning. In Sec. IIVI we give the applications of 
our SU(2) Kubo formula to some models of spin- 1/2 rep- 
resentation. In Sec. our theory is applied to a high- 
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dimensional representation (i.e., spin-3/2 representation) 
system. Several concrete example models are also given. 
In Sec. IVII we give a brief summary with some remarks. 
In the appendixes, we give the detailed calculations of 
the correlation functions in Matsubara formalism. 



II. SU(2) KUBO FORMULA AT ZERO 
TEMPERATURE 



A generalized total Hamiltonian for a variety of models 
to study spin transport can be written as [18j 

H = 2^ E(P' " l A (ri,t) - vA a (r h t)r a y+ V dis , (1) 



where Vdis is the potential caused by disorders. Through- 
out this paper, the index I refers to the Zth particle, a, b, 
and c refer to spin space while i and j the spatial space, 
and repeated indices are summed over. r a stands for 
the generators of SU(2) algebra, A and A a are the U(l) 
and SU(2) gauge potentials, respectively. Usually these 
gauge potentials consist of two parts, internal and exter- 
nal fields. In order to derive a general formula for the 
conductivity in response to an SU(2) external "electric 
field ", we separate the Hamiltonian Q into two parts, 
H = Hr, + H' with 



H = 77-^2^ + Vdis, 

i 



(2) 



where the operator ni = pi - §A int (r ; ,i) - rjAf nt (ri, t)r a 
stands for the dynamical momentum involving internal 
U(l) and SU(2) potentials if any. Note that this separa- 
tion is up to the first order of A£ xt (r/ , t) . 

The relation between SU(2) "electric fields " and the 
gauge potentials is given by 



Et = -d A? - diA% + r,e abc AlA°, 



(3) 



where r\ is the coupling constant. The SU(2) "electric 
field " Ell i s expected to be realized by spatially [2l| 
or timely [23 dependent Rashba or Dresselhaus coupling 
strength. We consider the linear response with respect 
to the "electric" components of a non-Abelian field, Ej = 
E°-T a . For simplicity, we take the SU(2) external field to 
be of single-frequency, namely, 



(4) 



The problem involving the external field of an arbi- 
trary form on r only differs from this case by a fourier 
transform. As is well-known, the external field of 
frequency(q, u>) actually reads Bf(q, u>) cos(q • r — Lut) 



which is real valued. The expression in Eq. Q is conve- 
nient for calculation. Thus in the subsequent formulas, 
only the real parts have the physical meaning. 

Now we choose a gauge that Aq = 0, which corre- 
sponds to the temporal gauge in the U(l) case, then the 
SU(2) external "electric field" comes into the Hamilto- 
nian through the vector potential Af, 

At(v,t) = -E?(r,t) = -^(q,c) e ^ r — *. (5) 

VjJ 10J 

Hereafter, we omit the subscript specifying the external 
field for simplicity. 

Based on the definition for single particle, one can de- 
fine the SU(2)-current operator for the unperturbed sys- 
tem 

J Q (r) = ^-^[{r a ,*/}5(r-r0+5(r-r i ){r a ,7r i }], (6) 

where the curl bracket denotes anticommutator and the 
velocity operator of the Ith particle determined by the 
Heisenberg equation of motion, ni/m — [r/, Ho]/(ih) is 
matrix-valued for the SU(2) case. In terms of this cur- 
rent, the perturbation Hamiltonian H' can be expressed 
as 



H> = - J dr J?{v)Al{v,t) 



Jf(q) £?(q,u;)e- 



(7) 



where Jf (q) is the Fourier image of Jf (r). Clearly, the 
interaction term is the product of the SU(2) current and 
the SU(2) "electric field 

Taking the perturbation of the external field into ac- 
count, we have 11/ = jt/ — r]A a (ri)T a . Then the total 
SU(2) current driven by the external SU(2) "electric field 
" reads 



—y 

i 



[K^/j^r-rO + ^r-rOK,^} 

n 2 

(8) 



= J a (r)-^ a (r,i) n , 

where no is the particle density. 

At zero temperature, this SU(2) current is evaluated 
for the ground state of the system. In the interaction 
representation, the state \ip(t)) of the system at time t 
is related to the eigenvector \<p) of Ho by the S-matrix, 
i.e., \ip(t)) = S(t,—oo)\<p). Up to the linear order of 
H'ji?), S(t, -oo) = 1-| dt'HW), where H'^t') = 
e iH t'/h H i e -iH t'/h_ Thcn the average of tiic nrst term 

in the total SU(2) current Eq. © is given by 



<Jf(r,i)) = (J«(r,t)) -1 J Jt'dJ^r^HW) ]) 
= I' dt'e^-^e-^H J?{r,t), J, fe (q,t')])o, 



(9) 
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where (jf(r,i))o has been dropped since no SU(2) cur- 
rent is considered to follow in the absence of the external 
fields. Together with the second term, we obtain the fol- 
lowing expression: 



(^(r,i)> = <^M)> 
E b Ar,t) rt 



Am 



A?(r,t) n 



= W ;r)£$(r,t). 



(10) 



Since the conductivity represents the property of the 
whole system, we need take the average over the system 
to get the SU(2) conductivity, 



1 



huv 



4muj 



dt> e -(*- t ')([jf(q,t),j}(q,<')]>0 

(11) 



with V the volume of the system. The spin conductivity 
here is a tensor in spin space rather than a vector as in 
the case of linear response to the U(l) external field. 

As a conventional strategy, a retarded current-current 
correlation function is thus introduced to calculate this 
conductivity, 



Q?,-(q,*-0 



V 



6(t-ti)([j! a ( q ,t), J*(q,i') ]) , 



(12) 

where 9(t — t') is the step function which vanishes unless 
t > t' . The Fourier transform of Eq. H12|l is given by 



Q#(q,w) 



V 



dt6{t-if)e i ^- i '\[j} a {c i ,t), JHq,t>) ])„. 



Comparing with Eq. Ijll|l . we obtain 



< b (q,c) 



fiuj 



4m 



S a0 5 l 



(13) 



(14) 



To simplify the calculations, we introduce a Matsubara 
function Q^-(q, iv) which reduces to the retarded corre- 
lation function Qfj(q,uj) by changing iv to u> + iS, 



III. AN EQUIVALENT FORMULATION FOR 
ZERO FREQUENCY 

In the previous section, we derived the SU(2) Kubo for- 
mula choosing the gauge potential Aq = 0. To obtain the 
dc conductivity, one just needs to take the limit ui — > 0. 
As is well-known in the conventional electrical conductiv- 
ity, the Kubo formula can also be derived alternately by 
choosing a constant external field as a start point. The 
continuity equation for electric charge conservation guar- 
antees the two derivations to be equivalent. Whereas, in 
the SU(2) case, the current defined by Eq. © is not con- 
served as long as an SU(2) interaction is present. For 
example, the spin current, a special SU(2) current with 
r] = H, is not conserved if there exists the Zeeman term 
or spin-orbit coupling. In this case, the continuity equa- 
tion does not hold 0,0], instead, we have the following 
relation: 

(—-r]Aox)a(r, t)+(—+r ] Aix)j i (r, t) = 0, (16) 



where a a (r) = T]^(r)T a tp(r) and J a (r,i) are the SU(2) 
density and current respectively, and notations a = 
{a 1 , cr 2 , a 3 ),A = (A 1 , A 2 , ^ 3 )etc. are adopted. Unlike 
the charge current which is conserved, the spin current 
is not conserved, thus a natural question is whether the 
SU(2) Kubo formula we derived in the previous section 
is still consistent with the other derivation? 

Now let us choose doAf — for the zero frequency 
case, then the SU(2) electric field and the perturbation 
Hamiltonian are given by 



Ef = -diAl + ye^A^A 



and 



H' 



dr cr a (r, t)A a {r) 



(17) 



(18) 



By means of the method suggested by Luttinger, the 
total SU(2) current can be obtained once the density ma- 
trix p is introduced. The density matrix including the 
deviations caused by the perturbation takes the form 



P^) = Po + 



(19) 



where p refers to the density matrix with respect to the 
unperturbed Hamiltonian and 5p(t) is brought about by 
the perturbation one, H' . From the equation of motion 
for the perturbed part of the density matrix, 



Q#(q>) 



V 



du e wu (T u 4 a (q, U )J, 6 (q,0)), 



(15) 

where T u denotes the it-ordering operator and (3 — 
(fc B T) _1 with k B the Boltzmann constant. We thus have 
derived a generalized Kubo formula for spin transport in 
response to an external SU(2) "electric field " . 



ih^- = [H ,Sp(t)} + [H',p ] 



we can obtain a solution for 8p(t) 



(20) 



i r°° r 13 d 

Sp(t) = -j: J o dt dp' Pog- t Hi(-t - 0). (21) 
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With the help of the density matrix, the SU(2) current 
can be then evaluated by taking the average 



(#(r,f))=tr(p(t)#(r)) 







= -- I dt j d/3'tv Po -H'^-t-0) J?(r 



(22) 



where the equilibrium part of the current tr(p j"(r)) is 
assumed to be zero. The derivative of H'j with respect 
to time t is calculated as 



d t H' I {-t) = J dr d t a a (r,-t)A a (r). 



(23) 



Using the "continuity-like" equation (|16fl and integration 
by parts, we have 

dM-t) = J dv(pe ahc (Al a c - A\J%) - diJ?)AZ 

= - J 'dr E?J?(v,-t), (24) 

where we did not write out the arguments in the first line 
for simplicity. Substituting it into Eq. I|22() . we obtain 

(#(r,*)> 
I r°° r& 



o Jo 



= - dt dp' dt'tx Po E J fc (r',-t-z/3')Jf(r) 



(25) 



Consequently, the dc SU(2) conductivity is obtained from 
the above equation after integrating r over the volume V, 



r ab 







d/3' tr p J b (-^-^/3') Jf . (26) 



This result is obviously independent on the frequency. 
It is also consistent with the one which we derived in the 
previous section once we introduce the representation of 
the eigenstates \n) of H . Note that the spin procession 
terms, r]AiXJi(r 7 t) — r/Ao X(?(r, t) , precisely compensate 
the second term of Eq. (|17l) , which makes our theory self- 
consistent. Since the SU(2) "electric field " includes an 
extra term of gauge potential in comparison to the U(l) 
field, the nonconservation of the SU(2) current plays an 
essential role in guaranteeing the consistency. That is 
to say, the SU(2) current exactly responds to the SU(2) 
"electric field " no matter which gauge is chosen. 



A. Spin susceptibility 

Spin is a category of SU(2) entity with ?; = h. The 
spin degree of freedom is discussed extensively in recent 
years for its promising application. The effective spin- 
orbit coupling emerged significantly in some semiconduc- 
tors [13)13 i s °f importance for its possible manipulating 
of spin. Using our SU(2) Kubo formula, we can directly 
calculate the spin susceptibility which describes the linear 
response of the spin density to the spin-orbit coupling. 

The spin susceptibility xf b is defined as 



(S a ) = X ?El 



(27) 



where S a = hJ2 k ClT a C k is the spin density. Here we 
adopted a simplified notion C\ = (CL, C|,) with 
creating a spin-up particle of momentum k etc. The cor- 
responding retarded correlation function in Matsubara 
formalism II° b (w) is given by 



n? 6 (if) = ~ 



'(T u S a (u)jt(0)). (28) 



Hereafter, we take the unperturbed Hamiltonian to be 
H = J2k C l( £ ( k ) + d a {k)r a )C k for its elegant form in 
Green's function. The second term represents the in- 
ternal SU(2) field with d a the components of this field. 
This system has two bands, E- = e(k) + \d\ and E+ = 
e{k) - \d\, with \d\ = ^d a d a . In the limit u -> 0, the 
susceptibility reads 



h \ - n F _ 



n F+ 



2V r 



Ml 3 



dk. 



(29) 



where np_ and n F+ are the Fermi distribution functions 
and "— , +" label the different bands. This result is anti- 
symmetric to the indices labeling spin degree of freedom. 
Using this result, we calculate the spin susceptibilities 
with two kinds of internal fields, Rashba and Dresselhaus 
couplings, 



H R = -2a(k xT y -k y r x ), 

h d = -2/3(k x r x -k yT y). 



(30) 



These two kinds of couplings dominate in narrow gap 
semiconductors such as GaAs and here we take their two- 
dimensional (2D) forms to represent the effective spin- 
orbit couplings in two-dimensional electron gas (2DEG) . 
In these cases, the components xl V vanish since d z = 0. 
The results are shown in Table I, where we have taken 
the usual parabolic form that e(k) — H 2 k 2 /2m. 



IV. APPLICATIONS FOR SPIN-1/2 
REPRESENTATION 



B. Spin conductivity 



From now on, we will give some applications of our 
SU(2) Kubo formula. In this section, we mainly focus on 
the examples in spin-1/2 representation without impuri- 
ties and in the limit q — > 0. 



With rj = H, the spin current reads 

J i ~2^ Ck \^k~ + dk~ T ' / 



(31) 
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SU(2) internal field 




\x 




A!/ 


Rashba (^rp—) 


1 








1 


Dresselhaus ( ) 





-1 


-1 






TABLE I: Spin susceptibilities: a and j3 are coupling con- 
stants for the Rashba and Dresselhaus coupling respectively. 



SU(2) internal field 




zy 

w XX 


21 

® xy 


/T 

^ xy 


21 

a yy 


21/ 

(T 

u yy 


Rashba (zrrfi — ) 


c 





c' 
2 


c 
2 





d 


Dresselhaus (-r^ — „) 





-c 


c 
2 


c' 
2 


-c' 






TABLE II: Spin conductivities: c and c' represent the dis- 
placements in k space. 



After calculating the Matsubara function (see Appendix 
A) and changing iv — > u>, we derive the conductivity 



ab 



1 \ ■> Up- 

2~V 



k 



\d\ 3 



n F+ c abc d c de ( k ) de ( k ) 

dki dkj 



(32) 



This expression manifests that the conductivity is anti- 
symmetric to the spin indices a, b and symmetric with 
respect to the spatial indices i,j with the parabolic dis- 
persion relation. Note that when the U(l) part of Hq 
is parabolic, i.e., e(k) = h 2 k 2 /2m, and d a is linear of 
k a , the summation over k vanishes. Since the conven- 
tional spin-orbit couplings are Rashba and Dresselhaus 
couplings, which contain no quadratic terms of fc, we 
should consider e{k) = ^ [(k x + c) 2 + (k y + c') 2 ] for non- 
vanishing results, which represents a shift of momentum 
k in the material. Table II shows the spin conductivities 
with two kinds of internal fields. 

At this stage, it is worthwhile to recall some previ- 
ous work in spin Hall effect. Up to now, a general con- 
sensus is made that the spin conductivity in response 
to an external Maxwell electric field is exactly canceled 
by the effect of disorder in two-dimensional electron gas 
with spin-orbit coupling. The cancellation is due to the 
parabolic form of unperturbed band structure |25j |. It 
is worthwhile to point out that our SU(2) conductivity 
also vanishes when e(k) takes the parabolic form even in 
the absence of disorder. It is an essential difference that 
our conductivity refers to the linear response to an ex- 
ternal Yang-Mills electric field which is also a vector in 
SU(2) Lie algebra space whose bases, the Pauli matrices, 
are anticommute. Anyway the system with nonparabolic 
dispersion relation is of great importance, since the con- 
ductivity, no matter in the usual spin Hall effect with 
disorder or derived by our SU(2) Kubo formula without 
disorder, is expected to be observed in experiments. 

Finally, we will consider an experimentally available 
case. Since the Rashba coupling strength can be tuned 
by the gate voltage applied to 2DEG, we take a Rashba 
coupling with time-dependent strength as the external 
SU(2) field and Dresselhaus coupling as an internal field. 



Then we can get an ac conductivity depending on the 
frequency w. The result reads 



XX 



M = - 



ch 6 o 



32to 2 /? 3 tt 



(33) 



where e D (uj) is the dielectric function caused by the Dres- 
selhaus spin-orbit coupling (2||, namely, 



4/3 3 



k 2 dk 



{2(3k) 2 - h 2 u> 2 



(34) 



This dielectric function is a macroscopic quantity and can 
be directly measured. Carrying out the integration over 
k gives a resonant result, 



T Z >J 



ch 2 



16nf3m 128tt/3 3 to 2 



In 



k-hw/2/3 



k + nw/2(i 



hp - 



(35) 



zx 

yy 



where hp— and refer to the Fermi momenta of both 
bands. The same resonance is also shown in Ref. p7l . 
Other components are given by a** — \<y z xy while a 
and cr^U differ from them by c — > c'. 



APPLICATIONS FOR SPIN-3/2 
REPRESENTATION 



In the previous section, we have discussed several ex- 
amples using the SU(2) Kubo formula in spin-1/2 rep- 
resentation. It is well- known that there exit many im- 
portant systems which carry out the spin-3/2 represen- 
tation of the SU(2) algebra, for example, the Luttinger 
model [23l | containing the intrinsic spin-orbit coupling, 
bilayer systems |24| taking into account of spin degree of 
freedom, etc. Thus it is worthwhile for us to extend our 
discussion to high-dimensional representations, such as 
spin-3/2 representation. The examples mentioned above 
are also discussed, which may be instructive for the ex- 
periments. 



A. General consideration 

The spin-3/2 representation of SU(2) generators read 
/ 





V3/2 





v^/2 

1 






1 

x/3/2 
V^/2 



/ - 


iV3/2 





^1 


i\/3/2 





—i 








i 





-n/3/2 


V 





i\/3/2 


/ 


/ 3/2 









1/2 












-1/2 


:) 


• 


\ 





-3/2 y 





(36) 
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To simplify the calculations, we use a convenient ret 
resentation of the Clifford algebra adopted in Ref. [28 



P 5 



namely, r 1 = o z <g> a v , T 2 = a z <g> 
r 4 = a x ® I and r 5 = a z <g> a z where the cr's are the 
2 by 2 Pauli matrices. These gamma matrices satisfy 
{T^P 3 } = 2(5 Q/3 and r 1 r 2 r 3 r 4 r 5 = -l. Hereafter, the 
Greek indices a,/3, 7, etc. run from 1 to 5. This rep- 
resentation can be obtained from the Dirac representa- 
tion of the gamma matrices by a unitary transformation. 
Those five matrices can compose ten antisymmetric ma- 
trices Y a @ — [r Q ,r ,9 ]/2i which constitute the spinor rep- 
resentation of SO (5) algebra. The spin-3/2 operators can 
be expressed as a linear combination of those T a P , i.e., 
r a = ^L?. fl r Q ' 3 , where the coefficients L a af3 = -L1 a are 



r a 



4i a/3 

just the five-dimensional representation of the SU(2) al- 
gebra (i.e., spin-2 representation of the angular momen- 
tum operators). 

Since T", r Q/3 together with the identity / span the 
space of 4 x 4 Hermitian matrices, one can write out a 
general Hamiltonian in spin-3/2 representation in terms 
of those gamma matrices, 



ff o =^^(e(fc)+d Q (fc)r Q )C fc , 



(37) 



where C| = {C], v C\ 2 , C| 3 , C| 4 ) with the second in- 
dex referring to either spin-band or spin-layer labels. 
Here we do not include the linear combination of T al3 
which makes the Green's functions difficult to calcu- 
late. For this unperturbated Hamiltonian, there exist 
two types of perturbation part H' . One is constructed 
by T a @ and the other by T a . The problem relating the 
spin current in Luttinger model in response to the effec- 
tive field caused by the structure inversion asymmetry 
is of the first type. In this case, the structure inversion 
asymmetry is taken as the perturbation and hence 



H' 



k 



(38) 



Then the linear response of the spin current to H' reads 



(J?) 



bub 



(39) 



In calculating the retarded correlation function Qf b (ii/), 
we will encounter 



tr(G(fc, iLu n )j?G(k, iuj n + iiy)r b ) 

-±tr(G(M<){^ + ^r a , L a ^}G(k,iu; n + iu) L^T""), 

I 



(40) 



where G(k, iw n ) is the Matsubara function for which the SU(2) Kubo formula is then 
detailed calculation is given in Appendix B. 

Since the traces of gamma matrices are always real, 
the appearance of double T a P matrices makes Eq. I|40f) 
real and the summation of the Matsubara function also 
gives no imaginary contribution after changing iv — > uj. 
This directly results in a vanishing spin conductivity. 

The second type of perturbation is constructed by P, 



Jf) = o ba h a . 



(42) 



The corresponding retarded correlation function is given 
by 



1 



H' = ^C\h a T a C kl 



(41) 







l6a M = -^/ du e ivu {T u Jp(u)f a (0)) . (43) 
v Jo 



After changing iv — > uj and taking the limit u> — > 0, we 
and later we will discuss some concrete examples. The obtain the dc conductivity 

I 



b a 



4hv 



1 



d$ L „A uf_ —riF, 



(44) 



B. Concrete examples 



the spin conductivity for a bilayer system undergoing the 



Now we are in the position to discuss two concrete 
examples with the second type of H' . First, we calculate 
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Rashba coupling along opposite directions, 

H = e(k) + aa z ® (k x a y - k y a x ) + & x ® I. (45) 

In this model the tunneling between the two layers is 
included in which £ accounts for the tunneling strength. 
We take the SU(2) flux to be the external perturbation, 
H' = 4> a z ® a z , i.e., /15 = cj). A direct calculation of the 
spin conductivity gives the following result: 



x5 



1 /V3 ft 2 (q 2 fc 2 + 2£ 2 ) \ fe ^- 

47ra v / a 2fc2 + £ 2 ^ 4 mQf2 



(46) 



In the limit f — > 0, of 5 reduces to — -r^-. 

As another example, we take the tunnelling term to 
be the perturbation, that is H' = ^'<j x ® I and Ho = 
e(k) + aa z ® (k x a v — k y a x ), correspondingly, h 4 = 
Then we have the following result: 



0. 



-X'4 



T .y4 



87ra' 



(47) 



VI. SUMMARY AND REMARKS 



In this paper, we have generalized the Kubo formula 
to describe the linear response of the SU(2) current to 
the external SU(2) "electric field", which traditionally 
describes the one to the U(l) external field. From two 
distinct routes with different gauge fixings, we derived 
the SU(2) Kubo formula and showed that those two 
approaches are equivalent. The non-Abclian feature of 
SU(2) electric field involves one more term of gauge po- 
tentials in comparison to the U(l) case, while this term 
precisely compensates the nonconservation part in the 
SU(2) continuity-like equation for the SU(2) current. 

For the interests in spin transport, we applied our for- 
mula to calculate the spin susceptibility and spin con- 
ductivity in the system containing a Rashba or Dressel- 
haus field. The results show that in the usual system, 
where e(k) — h 2 k 2 /2m, the spin susceptibility is con- 
stant. However, the spin conductivity vanishes, much like 
the case in the spin Hall effect where the spin conductiv- 
ity in response to the external electric field vanishes in 
the presence of disorder. To derive the nonvanishing spin 
conductivity, the systems with nonparabolic unperturbed 
band structure are necessary, and the spin conductivity, 
no matter in response to the U(l) or SU(2) electric field, 
is expected to be observed in such systems. What is 
more, we also discussed an experimentally available case. 
In response to the time-dependent Rashba field, the spin 
conductivity is related to the dielectric function which 
can be measured directly. 



Generalized to the high-dimensional representation, 
our SU(2) Kubo formula is available to discuss the Lut- 
tinger model as well as bilayer spin Hall effect. The spin 
conductivity in response to the effective field caused by 
structural inversion asymmetry in the Luttinger model 
always vanishes. 

The work was supported by NSFC Grant No. 10225419. 



APPENDIX A: SPIN CONDUCTIVITY IN 
SPIN-1/2 REPRESENTATION 

The correlation function of spin conductivity in Mat- 
subara formalism reads 

Qlhiv) = -I f du e^iTjf {u)J b M), (Al) 



VJ 



where 



tv-kZ^iT + m*-^ (A2 » 



After using Wick's theorem and introducing the Mat- 
subara function G(k, u) — — (T u Ck(u)Cl(0)) , one can ob- 
tain the correlation function 

k,LU n 



(A3) 

where G and G+ refer to G(k,iu>) and G(k,iuj + iv), re- 
spectively, which can be derived from the Fourier trans- 
form 



G(k, u) = ~ ^2 G{k, iu n )t 



In the case Ho = e(k) + d a r a , 
G(k,ioj n ) 



(A4) 



ihu> n + fx — Hq 
= f(k,iu n )(g{k,iuj n ) + d a T a ), (A5) 



with 



/(MWn) " (i^„ + M - E )2-|d|V4' 
g(k,iu> n ) — i?ki) n + n — e. (A6) 

In the last line of Eq. I|A5|I . G(k,iuj n ) is separated into 
two parts, the U(l) part and SU(2) part, which facilitates 
our calculation of the trace term, 
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tr 
2 

+ 



2~ ' Jl dh dk 
(±g(k, ioJ n )g(k, iu n + iv)5 ab + 2d a d b - d 2 



T , T 



<%,; Ski- 



dd- 2 
<9fc, 



/ \ <9d a <9d & i 

+ [g{k, iu> n )g{k, iw n + iv) + |d| 2 /4j — — j (A7) 



Note that summing the Matsubara function over the frequency gives 

~a X] -f( k > ^)f(k, iuJn + iu )9{k, iu n + if) = X] /(ki iuj n)f(k, iu n + w)g{k, ko n ) 

= \d\[(ihi,)*-\d\*y 



(A8) 



and 



~j3 ^2 iuJ n)f(k, iuj n + iv) = -T^p^y^ /(fc, iu n )f(k, iuj n + iv)g{k, iuj n )g{k, iiv n 



Ml 2 / 3 „.. 
2(riir_ — ?^f+) 
= |d|[(ife/) 2 -|d| 2 ]' 
Thus the last line of the right-hand side of Eq. (|A7(1 vanishes. Changing 

<5 ab |d| 2 - d a d b de{k) de{k) 



+ iu) 



(A9) 



parts of Qf b {u; 



2U^|dp 2 ^ 2 - |d| 2 ) fljfe^ 



ii> — > u>, we can obtain the real and imaginary 
(njr_- n F+ ), 



e abc d c 



de{k) de(k) 



2V ^ \d\(h 2 uj 2 - \d\ 2 ) dk t dkj 



(n F - - n F+ ). 



(A10) 



Since o$(w) 

the limit lj — 
Eq. Q2|, 



ft a; 



then taking also be expressed as the linear combinations of T a ^ , i.e., 



0, we can obtain the dc conductivity in 



Re a. 



a & _ 



lim — Im Q«*(u,) 



1 x rip- 
2U <^ 



nF+ abc , c de^de{k) 

\d\* 6 d ~dkT^k~- (A11) 



APPENDIX B: SPIN CONDUCTIVITY IN 
SPIN-3/2 REPRESENTATION 

Before we calculate the spin conductivity in spin-3/2 
representation, it is wise to warm up with the Clifford al- 
gebra. The 4 by 4 gamma matrices T a are constructed by 
the 2 by 2 sigma matrices, which satisfy {T a , T' 3 } = 25 a p 
and r 1 r 2 T 3 r 4 r 5 = — 1. Using these gamma matri- 
ces, one can also compose ten antisymmetric matrices 
pa/3 = ^Y a ,T^]. Together with the identity matrix, 
T Q and T a P span the space of 4 x 4 hermitian matrices. 
The SU(2) generators r a in spin-3/2 representations can 



r a = ±L a R T a P with 

4z ap 



L x 



V 



/ 








i 


iV3\ 







—i 










i 













-i 











\ 


-i-s/3 








o / 


( 





i 





o \ 










i - 


-iV3 




-i 













-i 
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iV3 








o / 
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-i 
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i 



















-2i 










2 1 








\ 











o / 



(Bl) 



Note that L a are antisymmetric and satisfy the commu- 
tation relation [L a ,L b ] — ie abc L c . Thus they form the 



spin-2 representation of the SU(2) algebra. The follow- 
ing formulas are inevitable in further calculations, 

tr(r a r /3 ) = 4<5 Q/3 , tr(r a r /3 r 7 ) = o, 

tr (r c T /3 r'T I/ ) = 4(S a/s S^- 5 a »5 l3v + S^S 13 ^), 

tr(r a r /3 r 7 r^r ,y ) = -a*?^ tr(r Q/3 ro = o, 

tr(T al} TT^T v ) = Aie^^. (B2) 

We take the unperturbed and perturbed parts of the 
Hamiltonian to be H = T,k C l( e ( k ) + rfa r Q )C fc and 
H' = Y<k C l( hl3TP ) C k- Accordingly, the Kubo formula 
for the spin conductivity reads 



(J b (v,t))=a ba (q,u;)h a (r,t) 
In the limit q — > 0, 

r t 

a uo(t-t') 



(B3) 



a b i a (u;) = j± v J Jt'e- '<[./('(/ ). F> (/')]}„. (134) 



and the corresponding retarded correlation function in 
Matsubara formalism is given by 



Q»» = -I [ P due^(Tj*(u)f a (0)) a . (B5) 



Similarly, introducing the Matsubara function G(k,iuj n ) 
and using the definition of spin current 



k 



we can calculate the trace term 



Ai 



-M 



(g(k, ico n ) + d"r") { ^ + ^r", (g(k, iu n + iv) + d^)Y a 

= 2{g(k, iu n + iv) g(k, iu> n j) ( 2 ^L^^ ~ e^i^^d") • 
Note that L a are all imaginary, then the dc conductivity is given by 

1 - 4 y^ im V 2 ^ af} ^dki) |d|3 ' 

I 



(B7) 



(B8) 
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